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Abstract 

We have developed a model to predict the post-collision brightness increase of sub-catastrophic 
collisions between asteroids and to evaluate the likelihood of a survey detecting these events. 


It is based on the cratering scaling laws of Holsapple and Housen (2007) and models the 


ejecta expansion following an impact as occurring in discrete shells each with their own 
velocity. We estimate the magnitude change between a series of target/impactor pairs, as¬ 
suming it is given by the increase in reflecting surface area within a photometric aperture 
due to the resulting ejecta. As expected the photometric signal increases with impactor 
size, but we find also that the photometric signature decreases rapidly as the target aster¬ 
oid diameter increases, due to gravitational fallback. We have used the model results to 
make an estimate of the impactor diameter for the (596) Scheila collision of D = 49 — 65m 
depending on the impactor taxonomy, which is broadly consistent with previous estimates. 
We varied both the strength regime (highly porous and sand/cohesive soil) and the tax¬ 
onomic type (S-, C- and D-type) to examine the effect on the magnitude change, finding 
that it is significant at early stages but has only a small effect on the overall lifetime of 
the photometric signal. Combining the results of this model with the collision frequency 


estimates of Bottke et al. (2005), we find that low-cadence surveys of ~one visit per luna¬ 


tion will be insensitive to impacts on asteroids with D < 20km if relying on photometric 
detections. 
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1. Introduction 


The main asteroid belt is collisionally dominated with large asteroids’ shapes, sizes 
and surface geology controlled by impacts. Studies of collisions help us to understand 
the evolution of the shape of the asteroid population and in turn the formation of our 
Solar system. These studies may involve laboratory experiments, computer modelling or 
observational programmes. 


The evidence for collisions can be seen indirectly in main-belt asteroid families (Cellino 


et ah, 2002), asteroid satellites and binaries (Merline et ah, 2002). It can also be seen 


directly in recently observed collisions (jSnodgrass et ah 2010; Jewitt et ah| 2011; Stevenson 


et ah, 2012). There are three possible collisions observed to date. In 2009 the 120 m 


diameter asteroid P/2010 A2 suffered a collision with a 6-9 m estimated diameter impactor 
( Snodgrass et ah[|2010 ) (but see section 4.4). In 2010 another asteroid, (596) Scheila (113 
km diameter), was hit with a ~ 35m diameter impactor Jewitt et ah (2011). The most 
recent potential collision involved the object P/2012 F5 (Gibbs), which like others was 
originally identified as a potential main-belt comet Stevenson et ah (2012). Events like 
the (596) Scheila collision should occur approximately every 5 years and collisions with 
asteroids <10m even more often (Bodewits et ah, 2011). 

Several recent surveys are capable of detecting collisions and cratering events. For 
example, the Canada-France-Hawaii Telescope Legacy Survey was used to search for Main- 


Belt comets among 25240 objects in 2003-2009 (Gilbert and Wiegert, 2010), the Thousand 
Asteroid Lightcurve Survey (924 objects) was conducted with the Canada-France-Hawaii 
Telescope in September 2006 (Masiero et ah, 2009) and the Hawaii Trails project was 


conducted in 2009 (599 objects) (Hsieh, 2009). While none of the surveys mentioned above 
were specifically looking for main belt collisions, the methods used in search for main belt 
comets would have also revealed any collisional events. There are also current surveys 
fully or partly dedicated to discovering Near Earth Asteroids, such as Pan-STARRS 1 
( Kaiser et ah| 2002), the Lincoln Near-Earth Asteroid Research (LINEAR, responsible 
for discovery of P/2010 A2) (Stokes et ah, 2000), the Catalina Sky Survey ( |Spahr et ah 


and the VST ATLAS survey (Shanks et ah, 2013) that are all capable of detecting 
main-belt collisions. 

Much work has been done in modelling the parameters (i.e. shape of debris, brightness. 


total ejected mass, impactor mass) of known observed collisions (Kleyna et ah, 2013), 
(Ishiguro et ah, 2011a), (Holsapple and Housen, 2007), (Housen and Holsapple, 2011) ; 
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and hydrodynamic modelling of generalised collision (Benz and Asphaug, 1999). This 


work focuses solely on the magnitude change following an impact as it is most likely to be 
observable by optical telescopes. Rather than looking at a specific object in the main belt, 
the described model looks at what would be expected with generic asteroids. 


2. Model description 

2.1. Cratering physics 


Our model is based on the work by Holsapple and Housen (2007), who provide a 


summary of scaling laws that allows calculation of crater size using properties of the target 
and impactor, based on the results of impact experiments. These laws can also be used 
to calculate the evolution of the ejecta dispersal and consequently estimate the amount of 
material ejected and increase in brightness following a collision. The decrease in magnitude 
of the target asteroid is going to depend on the amount of material that was ejected and 
whether it is optically thin or not. 

At high impact speeds, transfer of the energy and momentum of the impactor into the 
target occurs over area on the order of impactor size, while the resulting crater usually 
exceeds this size by many times. It is therefore a reasonable approximation to assume that 
impact occurs as a point source. Using theoretical analyses of mechanics of crater forma¬ 
tion, Holsapple and Housen showed that the crater and ejected material characteristics 
depend on the quantity aU^5'^, where a is the radius, U is the normal velocity component 
of the impactor and <5 is the density of the impactor; fi and v are scaling exponents. 

The scaling exponents depend on the material properties. Theoretical values of p, range 
from 1/3 to 2/3 ( Holsapple and Schmidt| 1987) and are a measure of the energy dissipation 
by material; a more porous material can dissipate energy more effectively and will have a 
lower value of this exponent. Experimentally determined values of p are ~0.55 for non- 
porous materials (e.g. rocks and wet soils), 0.41 for moderately porous materials (e.g. sand 


and cohesive soils) and 0.33 to 0.40 for highly porous materials (Holsapple and Schmidt 


1987). Experimental values for u were found to be the same for all materials at around 0.4 


(Holsapple and Schmidt, 1987). By selecting appropriate material scaling parameters for 


a given impact and inserting them into a general expression for the relationship between 
radii of involved objects and crater size, a reasonably accurate estimate of the crater size 
(as well as crater formation time and transient crater growth) can be made. 


3 















We now summarise how we use the previous studies in our calculations. Consider a 
spherical, non-rotating asteroid of radius r following an impact from an object with radius 
a at sub earth point. The general form of equation for crater size R consists of strength 
and gravity term: 


R = aKi (gravity term -|- strength term) 


. A* 
2 +^ 


( 1 ) 


where 


So-a, p . 21 ; 
gravity term = -) '' 

Op- 


•^grain 


strength term = 


Y , 2 +^ , 8, 






( 2 ) 


( 3 ) 


>U2' ■ p 

Here Ki is a scaling parameter (1.03, 1.17, 0.725 for sand/cohesive soil, wet soils/rock 


and highly porous material respectively; (Holsapple and Housen, 2007)), Y is the average 
strength of the target material, p is the grain density of target, U = 5 km s“^ (Bottke 


et al., 1994) is the normal velocity component, 6grain is the grain density of the impactor 


Sg is the surface gravity of the target asteroid with mass M and radius r, calculated as 
follows: 


s,= 


GM 


( 4 ) 


Depending on the asteroid type, different values of bulk density (for calculation of target 
asteroid mass) and grain density (for calculation of ejected mass) are used. The values and 
their sources are summarised in Table Bulk densities of C- and S-type asteroids were 
taken from weighted averages of corresponding subclasses as summarised in Table 3 of 


Carry (2012). Grain densities of C- and S-types are assumed to be the same as their most 


likely meteorite analogues (Britt et ah, 2002). Density of D-type asteroids is approximated 


by bulk and grain densities of the Tagish lake meteorite (Zolensky et al., 2002 Izawa et al. 


2010 ). 


The range of material strengths used is presented in Table The strength value selected 
in this study was varied for each taxonomic type to explore the relationship between the 
type, strength and the corresponding magnitude change. 

The crater radius R calculated in this way has a corresponding mass Mcratei- 
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Table 1: Average bulk and grain density of asteroids depending on taxonomic type. 
Asteroid type Bulk density (kg m“^) Grain density (kg m“^). 


C 

s 

D 


1840 ( 

Carry 

201^ 

!) 

2640 ( 

Carry 

201^ 

!) 

1670 ( 

Zolensky et al. 

2002) 


2710 

Britt et al. 

2002) 

3700 

Britt et al. 

2002) 

2770 ( 

zawa et al. 

2010) 


Table 2: Average material strength used in calculations for a selected scaling. 


Material Average strength Y (MPa) 


Sand/cohesive soil 

Wet soil/rock 

0.18 ( 
1.14 ( 

Holsapplel 

Holsapplel 

1993) 

1993) 


Highly porous 0.001 ( 

Holsapple and H 

ousen 

2007 

) 


M, 


crater 


= k. 


crater 


pR^ 


(5) 


where the scaling factor fecrater is taken to be 0.75 for cohesive soils, 0.8 for wet 
soils/rocks or 0.4 for highly porous material (Housen and Holsapple, 2011j ). 

As we are interested in the ejected mass, since it is only that which contributes to the 
observed magnitude change of the asteroid, the full crater mass will give an overestimate 
of brightness. The total crater volume is made up of a volume of ejected mass, a volume 
of the mass that is uplifted near the crater rim and a volume due to compaction. The 
fraction /^ejecta of the total crater mass that corresponds to ejected mass is of order 0.2-0.5 


(Housen and Holsapple, 2011) 


Afejecta — ^ejectaAfcrater (6) 

Throughout this work we assume fcejecta of 0.3 as being most appropriate to asteroids. 
2.2. Velocity shell model 

We consider the ejecta leaving the asteroid surface after the collision event. For sim¬ 
plicity, we assume that the debris expands spherically outwards from asteroid, the debris 
with each velocity forming a shell of radius (see Figure [^. Effects from rotation 
of the target or impactor are beyond the scope of the current model. Impact experiments 
show that there is no significant correlation between velocity and mass of the particles 
(Holsapple et al. 2002). Therefore, each velocity shell is taken to have the same parti- 
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cle size distribution described below, in Section 2.3 As our aim is to model observable 
brightening from Earth, we assume that the ejecta cloud is centred on the asteroid, as at 
early epochs the asteroid itself and the ejecta will be unresolved. We also assume that 
the brightness of the asteroid plus ejecta is measured through an aperture of fixed radius 
Tap and centred on the asteroid. To obtain the amount of material that is visible in the 
aperture, and consequently the visible increase in asteroid magnitude, we need to look at 
the fraction of each debris shell that fits completely or partially in the aperture (ejecta 
visibility fraction /vis)- 



(a) Velocity shells. 



Figure 1: Velocity shell model: (a) Ejecta expands spherically outwards from the asteroid. The shell 
expands with velocity u„. The portion of each shell that fits into the aperture rap is marked in solid line, 
(b) Spherical cap indicates the portion of each shell that fits the aperture. 


There are three possibilities for an individual shell at a given time t: 
1. the entire shell fits within the aperture and hence /vis = 1 
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2 . none of the shell is visible, i.e. the material is not above the surface of the asteroid 
as it has fallen back, /vis = 0 

3. part of the shell fits into the aperture (0 < /vis < 1) (see Figure 

In the latter case, the amount of material visible can be estimated by calculating the 
ratio of the surface area corresponding to the arc that fits into the aperture (spherical cap) 
to the total surface area of the sphere. 


/vis = 


2{2'iTrh) 

471 - 7-2 


( 7 ) 


Here h is the height of the cap and the extra factor of 2 accounts for near and far side 
material, h = r — ^ 7-2 _ 7 . 2 ^ 7 -^^ jg aperture radius that equals the radius of the 

cap. In the situation of no gravity acting on the debris, the radii at each velocity can be 
described straightforwardly by Vg = vt, however we need to consider the gravitational case. 
For simplicity, we consider a single debris shell of radius and velocity v over the asteroid 
of radius r. The acceleration due to the asteroid’s gravity that the debris experience is 
inversely proportional to the square of shell radius (-^) and also happens to be the second 
derivative of radius with respect to time. 


(frs 


Cr 


-2 


( 8 ) 


C is a constant of proportionality which we calculate by considering the boundary 
conditions. At the asteroid surface = r and the acceleration due to gravity is the surface 
gravity Sg. At the surface the original equation reduces to Cr~^ = Sg, which in turn gives 
C = Sgr‘^. Therefore, the radius of the spherical shell of debris at velocity v is governed 
by a second order differential equation of the form: 




dt"^ 



( 9 ) 


with the boundary conditions that at t = 0, rg = r, -Jf = ^initial- Solving this equation 
will give the shell radii at all times t. The total mass within the aperture Mvis is the integral 
of all velocity shell fractions within the aperture i.e. the shell mass Mgheii multiplied by 

/vis 


Mv 



Alshell/visdu 


In the code it is calculated by taking a sum of the individual contributions from each 
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velocity shell, thus 


A^vis — 5]^™^^.^Mshell/vis 

As the asteroid fragments will be ejected at a continuous range of velocities, increasing the 
number of shells considered improves the accuracy of calculations. 


2.3. Magnitude change calculation 

At this stage we implicitly assume that the scattering function for the asteroid and the 
ejecta particles is the same and that the ejecta is optically thin throughout. Assuming all 
the ejected particles contribute to reflected light with the same albedo as the target asteroid, 
the total visible surface area after the collision is Ac = Aa + Ap where Aa is the area of 
target asteroid, Ap is the area of ejected particles. The latter can be calculated assuming 


that the particle size distribution follows a power law of the following form (Ishiguro et al. 


2011 a): 


J.J/ \ ] f ■^o(o„) da, 0,min ^ ^ 0,max 

N{a)da = < ° 

I 0, O < Clmim 0. > 0-max 


( 10 ) 


Here Omin = 0.1 /rm is chosen to be the minimum particle radius because scattering is 
inefficient for particle sizes much smaller than the wavelength of the scattered light, Omax 
is the maximum particle radius, Nq is the reference dust abundance at reference size oq 


and q = 3.5 is the assumed power law index (Jewitt, 2009). The maximum particle radius 


Omax is calculated using the expression derived from the power-law function of the ejection 


velocity of dust particles (Ishiguro et ah, 2011a) 


Ojm.n.x — 



( 11 ) 


Here Vo=80 m s“^ is the reference ejection velocity of 1 /rm particles, k = 1/4 is the 
power index of size dependence of ejection velocity, M and r are mass and radius of the 
target respectively. 

After integration we get cross-sectional area per unit mass of Oa (m^ kg“^), and the 
total cross-sectional area of particles around the asteroid Ap are: 


,-0.5 


-0.5 


Oa = - 


4p, 


^ ^max 

0.5 


■gram ^max 


aO-5 - a 


( 12 ) 














Ap — (Toiikfvis (1^) 

The change in magnitude is going to be related to the ratio of the total area after 
collision Ac and initial target asteroid area Aa- 

Am =-2.51ogio ^ (14) 

■^a 

Any ejecta in front of the asteroid will obscure a cross-sectional area equal to its own 
area. Given that we assume the same scattering function for asteroid and ejecta, this means 
ejecta directly in front of the asteroid does not contribute to the net brightness. Any ejecta 
behind the asteroid is clearly not visible. Therefore all ejected material along the line of 
sight of the asteroid should be ignored for brightness calculations. This is implemented by 
the subtraction of a pair of spherical caps with radius equal to that of the asteroid. 


3. Modelling results 


3.1. Model setup 

We have looked at three taxonomic types that are common in the main belt - C, S 


and D types (DeMeo and Carry, 2014). The differences in the types are accounted for 


using their corresponding grain and bulk densities during calculation of the ejected mass 
and magnitude change. In the simulation we varied the taxonomic type (C, S, D) and 
strength (sand/cohesive soil and highly porous) value to see if any patterns emerge. We 
have used a 50 by 50 point grid of target and impactor radii, ranging from 1-100 km and 
1-100 m respectively and logarithmically spaced, and 500 velocity shells linearly spaced 
over interval 10“^ to 10^ m s“^. The lowest initial velocity has to be above zero for the 
purposes of calculations and the amount of material being ejected at and beyond 10^ m 
s“^ is a very small fraction of the total ejected mass and therefore considered negligible, so 
we take 10^ m s“^ as a upper limit of the initial velocity range. Given that there is a range 
of target and impactor pairs that would result in catastrophic disruption (i.e. a dispersive 
shattering event leaving no remnant larger than 1 /2 of the original mass of target asteroid 


(Greenberg et ah, 1978)) rather than in sub catastrophic collision, the results of our model 


for that region will not be physical. We have estimated the location of the disruption region 
by calculating the approximate impactor radius necessary for catastrophic disruption of a 
basalt target assuming average density of 2500 kg m“^ and impact velocity of 5 km s“^ 
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(Benz and Asphaug, 1999). This region is marked in black on the figures. The code for this 


model was written in Python 2.7 (includes NumPy and SciPy packages), plotting of the 
output was performed in MatLab. The code requires a separate run for each pair of values 
for the taxonomic type and strength regime (i.e. C-type in sand/cohesive soil and C-type 
in highly porous regime would be two different code runs). Each run from impact until 7 
days after impact with the above parameters takes approximately 2 hours on a single core 
of an Intel Core i5 2410M processor running at 2.3 GHz. 


3.2. Model Results 

Figure shows the resulting predicted change in magnitude post-collision for a C-type 
asteroid with both a highly porous and a sand/cohesive soil strength at two epochs after 
the impact (1 and 7 days). Figureshows the calculated time it takes for the brightness 
increase to reach —1.0 magnitudes for C-types in the highly porous and the sand/cohesive 
soil strength regimes. This value has been chosen because the one sub-catastrophic aster¬ 
oid collision confirmed so far between the known main-belt asteroid (596) Scheila and an 


impactor was observed to have a brightness increase of ~ — 1 magnitude (Jewitt et al. 


2011). Conclusive attribution of a collisional nature to a smaller magnitude decrease may 


be difficult, due to natural variations in brightness caused by rotational modulation of the 


light curves and uncertainty in asteroidal absolute magnitudes (Pravec et ah, 2012). 
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Magnitude change 7 days post collision 
(C tai'get, C impactor, highly porous strength) 
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(a) highly porous regime, 1 day 


(b) highly porous regime, 7 days 
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Magnitude change 7 days post collision 
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(c) sand/cohesive soil regime, 1 day 


(d) sand/cohesive soil regime, 7 days 


Figure 2: C-type impactor collision with C-type target: 2(a)|2(by highly porous strength regime: magni¬ 
tude change 1 and 7 days post-collision for a range of target and impactor radii. 2(c)||2(d)| sand/cohesive 
soil strength regime: magnitude change 1 and 7 days post-collision for a range of target and impactor 
radii. Impactor radii range from 1 to 100 m, target radii: 1-100 km. Colour bar shows the magnitude 
change. Catastrophic disruption region is marked in black. Region where optical thickness of the debris is 
potentially significant is above the dashed line. Solid black lines indicate contours where magnitude change 
is -1, -3, -5, -7 and -9. 
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(a) highly porous regime, 1 day 


Time (days) at which magnitude difference 
falls to -1, truncated to 1 day 


(b) highly porous regime, 7 days 
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(c) sand/cohesive soil regime, 1 day 


(d) sand/cohesive soil regime, 7 days 


Figure 3: Time taken for the magnitude change to fall to -1 magnitudes after a C-type impactor collision 
with C-type target. Figures [3(a)| and |3(c)| show the impactor-target regime where this occurs within 1 day 
for highly porous and sand/cohesive soil strength regimes respectively. Figures 3(b)|3(dy show the same 
data but now extended to show where this occurs within 7 days of impact. For each figure the colour coding 
for this timescale is shown on the right, with the truncated region (timescale greater than 1 day or 7 days) 
shown in dark blue on the plot. The region where a catastrophic disruption would occur is marked in black. 
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Magnitude change after collision of a 20 m radius 
impactor with 1325, 1600 and 1930 m radius targets 



Figure 4: Magnitude change with time for three targets of radii 1325, 1600 and 1930 m after collision with 
20 m radius impactor (all S-type, sand/cohesive soil regime). Steps are clearly visible due to quantised 
nature of the modelled velocity shells and occur at a different position for each target. 
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There are vertical structures in these plots that are not of physical significance, but are 
rather features of the model limitations. Three processes affect the total magnitude de¬ 
crease post collision: the constant expansion of debris moving out of the aperture, material 
falling back on the surface of the target asteroid and therefore disappearing from observa¬ 
tion; and material re-entering the aperture due to the gravitational attraction of the target 
asteroid. In an ideal simulation of the process, the change in magnitude as a function of 
time would be a smooth function, being the result of the interplay of all three processes. 
The function would be smooth because we expect the size and velocity distributions of the 
debris to be continuous. As velocity shells in this model are quantised, this leads to the 
creation of artefacts at certain values of target radius. 

The effect of this can be seen clearly in Figure]^ which shows the predicted brightness 
increase after an impact of a 20 m radius impactor onto 1325, 1600, 1930 m radius target 
(all S-type, sand/cohesive soil strength regime). For lower velocity shells that fall back onto 
the asteroid at early times these steps merge to give a continuous slope. Higher velocity 
shells never return as they have reached escape velocity, and thus never produce a sharp 
decline. 

For a small number of shells with intermediate velocity, the return of the shell to the 
surface will be sufficiently distinct temporally from the other shells to produce a clearly 
identifiable step. The velocity shells leaving the aperture do not produce a sharp decline, 
because the fraction of the shell in the aperture decreases as a continuous function of time 
after the impact. The location of these steps is independent of the impactor radius and 
depends solely on the target radius due to the dependance on the surface gravity. These 
steps are clearly visible in Figures and after as vertical structures that should not be 
interpreted as being significant. 

One of the ways to minimise this and have a good resolution in the images is by using 
a large number of shells, however this increases the computation time. The 500 shells used 
in this study was selected as a compromise between resolution and run time. Figure 
shows a comparison between a plot using 500 shells and 5000 shells (20 hours run-time). 
The latter shows a reduction in vertical structures without substantial difference in output. 
Difference images show a discrepancy between the 500 shells and 5000 shells results being 
confined to the vertical structures only. 

We have chosen to only run the model up to 7 days post-collision, because a large 
proportion of the material (by surface area) would leave the aperture due to radiation 
pressure effects, which are not currently included in the model. The time t it takes for 
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(a) 7 days, 500 shells 


(b) 7 days, 5000 shells 


Figure 5: C-type impactor collision with C-type target, sand/cohesive soil strength regime: magnitude 
change 7 days post-collision for a range of target and impactor radii. Impactor radii range from 1 to 100 m, 
target radii: 1-100 km. Colour bar shows the magnitude change. Catastrophic disruption region is marked 
in black. Comparison of using (a) 500 and (b) 5000 velocity shells. The plot with 5000 shells shows a clear 
reduction in the vertical structures in comparison to 500 shells plot. Solid black lines indicate contours 
where magnitude change is -1, -3, -5, -7 and -9. 


material to escape the aperture radius rap is approximately equal to 



(15) 


where a is acceleration, that can be defined in terms of ratio between radiation pressure 
acceleration and the local gravity, /3, and gravitation acceleration to the Sun at 1 AU qq'. 
a = (^Qq. In turn, qq ss , where Rh is heliocentric distance. 

This gives the approximate time for radiation pressure to remove material from the 
aperture as 


(2rapRh^Y 

V I^GMq ) 


(16) 


where Rh is distance from the sun to the asteroid in question. 

Assuming (5 is on the order of ~ 0.1(Fulle, 2004) for the small grains in the ejecta 
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expected to dominate the scattering, ^ AU (mid-main belt), the time t is approxi¬ 

mately 2 days. Therefore, we do not present any model results that go beyond 7 days post 
collision, as without the inclusion of the radiation pressure they do not reflect physical 
reality. 


3.3. Optical depth effects 

Throughout the modelling, all ejecta within the aperture and not behind the asteroid is 
treated as contributing to the brightness increase via the projected surface area. Effectively, 
it is assumed that optical depth r is sufficiently low that the approximation r ~ 1 — e~'^ 
is valid. (For r < 0.2 this leads to a brightening overestimate of less than 10%). To 
estimate the potential optical depth effects, the ejecta cloud was divided into a series of 
concentric rings and the column density of each was calculated, then converted, using the 
particle size distribution, into the optical depth r. Any rings with r < 0.02 were ignored 
(this corresponds to an error of less than 1%), and any cases where more than 10% of 
ejecta mass was contained in rings with r > 0.02 were considered to have their brightness 
potentially overestimated. We found that at one day after impact, ignoring opacity may 
lead to overestimating the brightness increase by up to twenty percent (~ 0.2) magnitudes. 
The region where this may be relevant is outlined in shown in figures 2(a), 2(c)[ 7(a) and 
8 (a) I showing the brightness 1 day post collision. Clearly for all impacts apart from those 


on the largest asteroids, the correction to the calculated magnitude change will have little 
effect on detection and can be ignored. By one week after impact, we calculate that optical 
depth effects are negligible for all target impactor pairs. 


4. Discussion 

A collision between asteroids is a multi parameter problem requiring us to make some 
initial assumptions. This section is therefore divided into 4 parts. First we describe an 
application of the model to the only known collision in the parameter space explored - 
(596) Scheila. In the second part we examine the effect of the different strength regimes in 
a collision between C-type asteroids (as the most numerous in the outer belt where (596) 
Scheila is situated). In the third part we keep the strength regime the same (sand/cohesive 
soil) and vary the taxonomic types to examine the effect it has on the magnitude change. 
We conclude with application of the model to predicting the type of observable collisions 
by currently active surveys such as Pan-STARRS 1 and the Catalina Sky Survey. 
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4-1. Modelling the magnitude decrease in the (596) Scheila collision 

In December 2010 asteroid (596) Scheila was observed with the 0.68 meter Catalina 
Schmidt telescope to have increased in brightness by approximately 1.3 magnitudes in 


comparison to the previous month (Larson 2010). There have been multiple photometric 
observations of (596) Scheila reported as summarised in Table from which the impact 
is estimated to have occurred on December 3.5 UT, 2010 (Ishiguro et ah, 2011b). All 
observations measured the total magnitude, i.e. from the asteroid and the surrounding 


debris. Jewitt et al. (2011) used a 105,600 km aperture for the HST observations of 
(596) Scheila and this value is used for the calculations throughout this section. We have 
selected these observations for comparison with our model as giving the true magnitude 
decrease because they were taken by the same instrument and in the same observational 
circumstances. 


Table 3: Summary of observed change in magnitude for (596) Scheila 


Date of observation (UT) Approximate time Magnitude 

post-collision (days) change 


2010 Dec. 3.4 ( 

Larson 

201 

.0) 



0 

-1.3 

2010 Dec. 11.44-11.47 (Larson 

2010 

) 

7.9 

-1.1 

2010 Dec. 12 ( 

Hsieh et al., 2 

012 



9 

-0.86 

2010 Dec. 14-15 ( 
2010 Dec. 27.9 

Bodewits et a 

1. |2014 

) 

11 

-0.66 

1 Jewitt et al. 

2011 ) 


24.4 

-1.26 

2011 Jan. 4.9 ( 

Jewitt et al., 
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A collision with a smaller asteroid is the most likely cause of this brightening and there 


have been a range of impactor diameters suggested in literature: 35 m (Jewitt et al., 2011), 


30-50 m (Ishiguro et ah, 2011a), 30-90 m (Moreno et al. 2011), <100 m (Bodewits et al. 


2011). The impactor diameters above are calculated using the estimated mass of high speed 
ejecta and the relationship m{v) ~ m* for v = 10“^Ui (Housen and Holsapple, 2011) where 
m{v) is the mass of ejecta above velocity v, mi is the impactor mass and Vi the impactor 
velocity. 

Scheila has semimajor axis a = 2.926 AU, eccentricity e = 0.164, and inclination i = 
14.7°, which places it in the outer main belt. The most numerous asteroids in that region 


and therefore the most likely impactors are C-type (DeMeo and Carry, 2014), although 


(596) Scheila itself is a D-type (Yang and Hsieh, 2011). The taxonomic type determines the 
densities used in the model which affect the final output parameter of change in magnitude. 
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Figure shows the post-collision magnitude change during 60 days for a collision between 
Scheila and C-, D- and S-type impactors of diameters reported in the literature. We also 


used the observations from Jewitt et al. (2011) to find the best fit size for each impactor type 


by matching the calculated magnitude change in our model to the observed lightcurve. The 
results are 59 m, 49 m and 65 m diameter impactors for S-, C- and D-types respectively. 
As has been mentioned in the previous section, we believe that in its current state the 
model gives best description of reality if we limit the results to 7 days post-collision. 
However, in the case of Scheila collision, which is an important real observed collisional 
event, we lack a sufficient number of data points of good quality in the first 7 days to 
make a meaningful comparison. However, the large aperture used in this study includes 
a large portion of the ejecta even after being significantly affected by radiation pressure, 
allowing comparison with our model. The model results are in approximate agreement with 
estimated impactor diameters by other authors. The differences most likely come from a 
variety of initial assumptions made by the different researchers, i.e. density, type of ejecta 
particle distribution, range of particle sizes and the estimated ejected mass. It is also clear 
from these results, that to the first order the taxonomic type of the impactor asteroid is 
not an important factor in diameter estimation. 

4-2. Effect of strength regime on post-collision magnitude change for C-type asteroids. 

We consider a collision where both impactor and the target are generic C-type asteroids. 
Figures 2(a) and 2(b)| show the magnitude decrease post-collision in highly porous regime, 
while Figures 2(c)| and |2(d)| show the effect of changing the regime to the sand/cohesive soil. 
Following the collision, a shock wave travels through the target asteroid and the outcome 
of the collision depends on how the stress wave propagates and is attenuated through the 
target. The cratering and amount of ejecta following an impact is determined by the target 
internal structure, porosity and strength. A porous asteroid would react to the impact 
differently than a more solid body and it is particularly important to consider impacts on 
these as most observable main belt asteroids should have significant macroporosities, apart 
from the very largest. It has been found that highly porous objects such as (253) Mathilde 
are particularly good at attenuating stress waves due to energy going into the compaction 
of the pores ( ]Britt et al. 2002). An object with large macroporosity is most likely to be 
highly fractured or even a rubble pile. Dark and primitive asteroids are more likely to 
have significant porosity, however it is worth noting that it is not exclusive to C-types. 
Figures [ 2 (a) [/|2(b')] and 2(c)[/2(d')]are similar for both strength regimes at each point of time 
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Magnitude change following a collision of (596) Scheila Magnitude change following a collision of (596) Scheila 

with C-type impactors of varying diameters. with S-type impactors of varying diameters. 



Time (days) Time (days) 


(a) C-type impactors 


(b) S-type impactors 


Magnitude change following a collision of (596) Scheila 
with D-type impactors of varying diameters. 



(c) D-type impactors 


Figure 6: Magnitude change as a function of time following a collision of C-, S- and D-type impactors 
with a range of diameters with asteroid (596) Scheila. Individual points show reported magnitude change 
by various authors. Dot-dash line shows our fit to Jewitt et al.’s data. The best fit estimates of impactor 
diameter with our model are 59, 49 and 65m for C-, S- and D-type impactors respectively. 


considered, however the sand/cohesive soil regime shows a larger magnitude decrease due 
to more material ejected, while in the highly porous regime some of the impact energy goes 
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into compaction. This tendency is also reflected in plots of amount of time taken for the 
magnitude change to decrease to 1.0 magnitudes above the pre-impact brightness (Figure 
[^. The brightness increase caused by the ejecta reaches the assumed observable limit 
of approximately 1 magnitude faster for the highly porous case than the corresponding 
impactor-target pair undergoing a collision in the sand/cohesive soil regime. This is due 
to less ejecta being produced. 

Recent research by Carry (2012) indicates that the density may be dependant on target 
size, rather than being a fixed quantity for each taxonomic type. However, the data is sparse 
in the size range used in this model and unavailable for the D-types. Therefore, our model 
makes a simplifying assumption that the bulk density stays fixed for all sizes considered 
and depends only on the taxonomic type. 


J^.3. Effect of taxonomic type (C, D and S) on post-collision magnitude change in sand/cohesive 
soil regime 

It is important to understand the taxonomy effects on the collision lightcurve due to 
the established composition gradient across the main belt. To investigate we kept the 
strength regime (sand/cohesive soil) constant and varied taxonomic type of both target 
and impactor (for S-, C- and D-types). Figures and show the magnitude change 1 
and 7 days post-collision for target-impactor pairs. Other combinations of target-impactor 
types are included in the Supplementary material (available online). S-types impacting 
each other are more likely to occur in the inner belt, while D-types impacting each other 
are more likely in the Trojan population. Generally, both sets of plots show qualitatively 
similar results for magnitude change with time, with a high magnitude change in the day 
1 , followed by the decline over time, with C-types having the least magnitude increase at 
7 days. From our results we observed that collisions between the lowest grain density and 
highest porosity C-types produce the least amount of magnitude decrease lasting for less 
time that it does for other taxonomic types. The S-types which have lower porosity and 
the highest considered grain density show the opposite - a larger magnitude change lasting 
longer. D-type collisions fall between these regimes. 

Figure shows a direct comparison between different combinations of target and im¬ 
pactor taxonomic types for a 20 km target and 50 m diameter impactor. In a collision with 
an S-type target, the largest magnitude change is produced by a S-type impactor, then 
C-type and D-type producing the least magnitude change. The same pattern is followed 
by a C-type target being impacted by S-, C- and D-type asteroids. However, a D-type 
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target deviates from this pattern, with an S-type impactor giving the largest change in 
magnitude, followed by D-types and C-types giving the least magnitude change. Overall, 
S-type impactors produce the brightest collision signature in all types of targets considered, 
due to their inherent high grain density. 



Target radius (km) 


Magnitude change 1 day post collision 
(S target, S impactor, sand/cohesive soil strength) 


Magnitude change 7 days post collision 
(S target, S impactor, sand/cohesive soil strength) 



lo'* io‘ ur 

Target radius (km) 


(a) 1 day 


(b) 7 days 


Figure 7: S-type impactor collision with S-type target, sand/cohesive soil strength regime: magnitude 
change 1 and 7 days post-collision for a range of target and impactor radii. Catastrophic disruption region 
is marked in black. Region where optical thickness of the debris is potentially significant is above the dashed 
line. Solid black lines indicate contours where magnitude change is -1, -3, -5, -7 and -9. 


Detection from current sky surveys 

The catastrophic collision rate in the main belt has been calculated from observationally 
constrained dynamical models as ~ 1 per year at a diameter D ~ 100 m (Durda et al., 1998 


O’Brien and Greenberg, 2003; Bottke et ah, 2005). The general hope within the scientihc 


community has been that such collisions would be detected via the on-going wide-held 
NEO surveys. However, several objects initially identihed as potential collisional disruption 


events are now suspected of being caused by rotational disruption due to YORP spinup 

e.g. P/2010 A2 (Agarwal et al. 

2013), P/2013 P5 (Jewitt et ah, 2013). Additionally, 

a recent study by |Denneau et al. 

(2014|) using 1.2 years of Pan-STARRS 1 data found 


only one plausible candidate of a collision event, and concluded that collisional disruptions 
of 100-m scale asteroids may be extremely rare. Hence observing much more frequent 
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Target radius (km) 


Magnitude change 1 day post collision 
(D target, D impactor, sand/cohesive soil strength) 





(a) 1 day 


(b) 7 days 


Figure 8: D-type impactor collision with D-type target, Sand/cohesive soil strength regime: magnitude 
change 1 and 7 days post-collision for a range of target and impactor radii. Catastrophic disruption region 
is marked in black. Region where optical thickness of the debris is potentially significant is above the dashed 
line. Solid black lines indicate contours where magnitude change is -1, -3, -5, -7 and -9. 


sub-catastrophic collisions may be a viable method for constraining the overall collision 
rate. 

To look at the likelihood of detecting these events, we first look at collisions with 100 km 


diameter targets similar to Scheila. According to Bottke et al. (2005), the size of impactor 
that would disrupt a, D = 100 km asteroid is D > 25km and such a disruption would 
happen every ~ 10^ years in the main belt. Using their CoDDEM model size distribution, 
a subcatastrophic collision impact by a H > 50 m impactor should occur once per ~ 25 
years, in rough agreement with the single such collision observed since the start of modern 
surveys in the 1990’s. The current Pan-STARRS and Catalina surveys are complete down 
to an apparent magnitude U ~ 20 (equivalent to ~ 2 km in the main-belt), but are 
sensitive to asteroids down to U ~ 22. The Bottke et al. disruption frequency at H ~ 2km 
is approximately one per 1000 years, requiring an impactor H ~ 60 m. Scaling by the 
CoDDEM model size distribution predicts a collision between a D ~ 2km asteroid and a 
D > 10m impactor should occur approximately once per year. As ~ 50% of the asteroid 
belt is visible at any time, current surveys should be able to detect a sub-catastrophic 
impact once every couple of years. 
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Magnitude decrease with time post-collision for different 
combinations of target and impactor taxonomic types. 



time (days) 


Figure 9: The effect of varying taxonomic type in a collision of 20km target and 50m impactor on magnitude 
change as a function of time. 


So why are these collisions between smaller bodies not found? In reality, observational 
surveys are affected by factors such as weather, moonlight and occasional losses due to de¬ 
tector gaps in the cameras. Perhaps most importantly, current surveys have a cadence for 
most individual asteroids measured in weeks rather than days. Using a detection aperture 
of ~ 1000 km equivalent to ~ 1 — 3 arcsec as appropriate to surveys such as Pan-STARRS, 
our model predicts a brightening of > — 1 magnitudes for D ~ 10m impacts on D ~ 2km 
for only 10-20 days. In reality, effects on the ejecta such as radiation pressure and Ke- 
plerian shear may shorten this timescale further. We can use the calculations presented 


by Denneau et al. (2014) (their figure 5) to estimate the likelihood of detection. With a 
maximum magnitude change of ~ 2 magnitudes and resulting rate of dimming of > 0.1 
magnitudes/day, the probability of the Pan-STARRS 1 survey detecting such a collision 
during the hrst 1.2 years of operation was < 0.05. Therefore we conclude that observational 
losses together with the rapid dispersal of the ejecta can signihcantly decrease the possi- 
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bility of photometric detection of sub-catastrophic impacts smaller that the (596) Scheila 
event, and this implies a current low probability of detection using automated photometry 
software such as MOPS (Denneau et ah, 2013). Finally, another additional factor may be 


inaccurate absolute magnitudes in current reference catalogues, as investigated by Pravec 


et ah (2012). They found that smaller asteroids were systematically fainter than expected 


from catalogues. Hence predicted magnitudes would be brighter than in reality, and any 
small brightness increase due to a collision could be masked by an incorrectly calculated 
residual magnitude. 

Of course, during the first few days to weeks the ejecta may be visible as an extended 
ejecta cloud around the target asteroid, and could be detected via direct manual observation 
or software algorithms designed for comet comae detection, as in the case of (596) Scheila. 
On the other hand, if the cadence of a survey is a week or less i.e. significantly less than 
the decay timescale within the aperture, the chance of photometrically detecting a small 


impact should be substantially higher. In Figure 10 we plot the predicted total F-band 
magnitudes for C-on-C and S-on-S collisions 1 day and 7 days after impact, for asteroids at 
opposition at = 3.5 AU. (Although S-types are predominantly found in the inner main 
belt, we calculate the predicted apparent magnitude at the likely maximum distance). 

First, it is clear that for the larger (but less frequent) impacts, the total magnitude will 
be relatively bright and would produce saturated images in large aperture surveys such as 
Pan-STARRS. However it may be possible to deal with this situation by either fitting to the 
wings of asteroid image point-spread function, or by recognising that an expected asteroid 
in the field was rejected in software processing due to its increased brightness. For fainter 
impact events, it is clear that Pan-STARRS with a limiting magnitude of F ~ 22 is able to 
detect such impacts throughout the asteroid belt, as long as the asteroid is observed soon 
after the collision. Importantly, the forthcoming ATLAS programme will have a nightly 
cadence over the visible sky and have a limiting magnitude of F ~ 20 ( |Tonry 2011). 
From comparing both Figure 10 with the brightness increases presented earlier, we find 
that ATLAS would be able to detect almost all of our studied collisions as well, but here 
the survey cadence should not be an issue. Therefore we conclude that current and future 
high-cadence all sky surveys should be able to detect many more asteroid collisions at early 
epochs. 
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(a) C-type target impactor, 1 day 


(b) C-type target impactor, 7 days 
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Figure 10: Visual magnitude following a collision. 10(a)|10(b) show magnitude 1 and 7 days post¬ 
collision for C-type impactor and target in the sand/cohesive soil strength regime for a range of target and 
impactor radii. 10(c)|10(d) show magnitude 1 and 7 days post-collision for S-type impactor and target in 
the sand/cohesive soil strength regime for a range of target and impactor radii. Impactor radii range from 1 
to 100 m, target radii: 1-100 km. Colour bar shows the visual magnitude. Catastrophic disruption region is 
marked in black. Region where optical thickness of the debris is potentially significant is above the dashed 
line. 
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5. Summary 


The model presented in this paper predicts the brightness increases caused by impacts of 


small asteroids on larger asteroids with radii > 1 km. We use the scaling laws of Holsapple 


and Housen (2007) and Housen and Holsapple (2011) to estimate the amount of ejected 


material. Our model separates the ejecta into discrete velocity shells, and assuming an 
ejecta size distribution calculates the magnitude change post-collision by calculating the 
surface cross-sectional area of the asteroid plus ejected material in a small photometric 
aperture. The scope of the model is limited to non-rotating asteroids whose ejecta particles 
do not interact with each other, however by extending the model results to the large 
apertures used for the reported ejects from the (596) Scheila collision, we find an estimate 
for the impactor size of 40 — 65 m (depending on taxonomic type) similar to previous 
studies in the literature. The model could be further improved by introducing particle 
light scattering laws, radiation pressure on the ejecta and the effect of rotation of the 
target asteroid. However this will lead to a significant increase the complexity of the 
model and will be developed in future. We believe our results are generic enough to be 
used to estimate the possibility of detection of such events by current automated surveys 
such as Pan-STARRSl and the Catalina Sky Survey. The model estimates that within 
the parameter space examined (impactor radius 1 — 100 m, target radius 1 — 100 km, 
sub-catastrophic collisions) , a magnitude change of less than —1 is observable by an 
automated survey like Pan-STARRS 1 with effective aperture radii of 1000 km for only 
10-20 days, which implies low probability of detection given the current low cadence for 
individual asteroids. However, detection may still be possible by direct manual observation 
or software algorithms designed for comet coma detection. 
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7. Supplimentary material 
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Figure 11: D-type impactor collision with S-type target, sand/cohesive soil regime: magnitude change 1 
and 7 days post-collision for a range of target and impactor radii. Impactor radii range from 1 to 100 m, 
target radii: 1-100 km. Colour bar shows the magnitude change. Catastrophic disruption region is marked 
in black. Region where optical thickness of the debris is potentially significant is above the dashed line. 
Solid black lines indicate contours where magnitude change is -1, -3, -5, -7 and -9. 
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Figure 12: C-type impactor collision with S-type target, sand/cohesive soil regime: magnitude change 1 
and 7 days post-collision for a range of target and impactor radii. Impactor radii range from 1 to 100 m, 
target radii: 1-100 km. Colour bar shows the magnitude change. Catastrophic disruption region is marked 
in black. Region where optical thickness of the debris is potentially significant is above the dashed line. 
Solid black lines indicate contours where magnitude change is -1, -3, -5, -7 and -9. 
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Figure 13: S-type impactor collision with C-type target, sand/cohesive soil regime: magnitude change 1 
and 7 days post-collision for a range of target and impactor radii. Impactor radii range from 1 to 100 m, 
target radii: 1-100 km. Colour bar shows the magnitude change. Catastrophic disruption region is marked 
in black. Region where optical thickness of the debris is potentially significant is above the dashed line. 
Solid black lines indicate contours where magnitude change is -1, -3, -5, -7 and -9. 
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Figure 14: D-type impactor collision with C-type target, sand/cohesive soil regime: magnitude change 1 
and 7 days post-collision for a range of target and impactor radii. Impactor radii range from 1 to 100 m, 
target radii: 1-100 km. Colour bar shows the magnitude change. Catastrophic disruption region is marked 
in black. Region where optical thickness of the debris is potentially significant is above the dashed line. 
Solid black lines indicate contours where magnitude change is -1, -3, -5, -7 and -9. 
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Figure 15; S-type impactor collision with D-type target, sand/cohesive soil regime: magnitude change 1 
and 7 days post-collision for a range of target and impactor radii. Impactor radii range from 1 to 100 m, 
target radii: 1-100 km. Colour bar shows the magnitude change. Catastrophic disruption region is marked 
in black. Region where optical thickness of the debris is potentially significant is above the dashed line. 
Solid black lines indicate contours where magnitude change is -1, -3, -5, -7 and -9. 
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Figure 16: C-type impactor collision with D-type target, sand/cohesive soil regime: magnitude change 1 
and 7 days post-collision for a range of target and impactor radii. Impactor radii range from 1 to 100 m, 
target radii: 1-100 km. Colour bar shows the magnitude change. Catastrophic disruption region is marked 
in black. Region where optical thickness of the debris is potentially significant is above the dashed line. 
Solid black lines indicate contours where magnitude change is -1, -3, -5, -7 and -9. 
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